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1 R-MODE OSCILLATIONS IN THE CONVECTIVE ZONE OF THE SDN 

1.1 INTRODUCTION. P-modes and g-modes have been actively studied 
since, it first appeared that the Sun might be oscillating in 
modes of global scale. But there is another type of global mode, 
called r-modes from their similarity to the Rossby waves in 
meteorology. These toroidal modes have very long periods 
comparable to the sun's rotation period, while the p- and g-modes 
have periods of minutes or hours, and their motion is controlled 
by the Coriolis force, whereas for the p- and g-modes the 
restoring forces are compression and buoyancy. 

Also the motion of r-modes is almost purely toroidal (non- 
radial) , while the p- and g-mode motions are mainly spheroidal 
(radial). Disturbance in the Sun would excite both toroidal and 
spheroidal components. We cannot continue to ignore the r-modes 
in solar physics. The r-modes are thought to be damped in the 
Sun. But damping does not remove them from observational 
interest. The solar envelope is in constant motion. Differential 
rotation drives eddies whose vertical component would be a source 
for r-modes. Lower angular harmonics (Fig. 1.1) would also be 
expected to be generated by the eddy motion of horizontal flows 
at the top of large convection cells and especially the cells of 
global scale. Any detection of r-modes would give an important 
new probe of the sun's interior rotation. 

1.2 EDDY VISCOSITY. We computed various properties of r-modes 
that would occupy the solar convection zone. Large scale motions 
in the solar envelope are non-adiabatic since they are immersed 
in a turbulent medium. 

If simple diffusion theory applies, an energy loss can be 
computed from the coefficient of eddy viscosity 77= 1/3 Q V con L, 
where V con = L/t c (mixing length/conv. overturning time). We 
calculated values of this using the solar convective zone model 
of Baker & Temesvary (1966). Following Landau & Lifshitz, viscous 
damping acting alone will change the energy of an r-mode in one 
oscillation period by 

e v - ) 2 

where the deformation tensor for large scale velocity is 



In spherical coordinates for 1=1 and both 1=2 cases, the 
summation of D^j terms give compact expressions which can be 
integrated over angle. Assuming that 77 depends only on radial 
distance, x=r/R 0 , 

E v = "2»6R 9 ^ydx n 6} (1.3) 
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where 5 is the term for the angular harmonic (Table 1.1). This 
approach requires the r-mode to have an oscillation period, P >> 
t c (coriv. overturning time), a long wavelength X >> L (mixing 
length), and a viscous energy dissipation E y << E (total energy). 

Special attention must be given to layers close to the base 
of the convective zone because of a possible anisotropic 
viscosity coefficient there. In the lower part of the convection 
zone, most of the energy would be lost to viscosity in a single 
oscillation period. Fortunately, this lower region is not very 
important to the mode as a whole; only about 10% of the undamped 
linear mode energy would lie below this region (at about 0.85 R Q , 
where the base of the CZ = 0.75 R Q ) . An estimate of the total 
viscous damping per period is one-third for the (1,1) mode. All 
other modes are damped more strongly. We can conclude the 
following if current models of the sun's envelope are realistic. 
When any transient disturbance excited r-modes in the CZ, the 
modes cannot retain their initial amplitudes for more than a few 
solar rotations without a fresh injection of energy at the proper 
phase (see Table 1.2 for damping rate) . 

One way to seek r-modes is to attempt direct detection of 
their horizontal motion. For this approach, it is best to go from 
a coordinate frame rotating with the sun to an inertial frame. In 
an inertial frame, the middle columns of Table 1.3 give the 
rotation rate and period of various modes, based on a convective 
zone rotating with a period of 25.7 days. Viewed from the frame 
rotating with the earth (synodic) at 365.3 days two rates may be 
of interest: the synodic rate and the repetition rate, at which 
phenomena caused by one r-mode repeats itself on the apparent 
solar disk. Table 1.3 shows the range of periods by listing the 
lowest and highest angular harmonics. The values for frequency 
and period apply to any surface feature varying linearly with 
displ ace-ment , such as velocity. Convection is likely to be 
modulated by r-modes. The overturning times of the large cells in 
the lower half of the convection zone are comparable to the 
oscillation periods of all r-modes with low 1 values. 

Also the length scales of these deep seated cells are 
similar to those of r-modes with moderate 1. Finally, convection 
has vorticity in the proper plane to couple with a toroidal 
velocity field. Thus the physical conditions for coupling are 
close to ideal. If convection is modulated on very large 
horizontal scales at time scales about one month, one should 
expect similar modulation in the various forms of solar activity 
like sunspots and flares. Wolff (1974) early suggested a relation 
between large scale convection and solar activity and since has 
provided additional observational material based on modes in the 
interior. Surface features would be affected more efficiently by 
the envelope r-modes, if they are indeed excited. This gives 
reason to search the solar observational record for evidence of 
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Fig. 1.1 Schematic fluid motion of r-modes with various values 
of 1 and m, near surface of the sun (after Saio, 1982). 


Table 1 . 1 Spherical harmonic terms of the stress tensdr for low 
order r-modes (after Wolff). r=radial dist; f=angular func(<^ in eg 1.3) 




* (l.i) 


■ -rfcos 0 


D «0 
rr 




r f coses i ne 



°00‘r(zrL dd 0 * V ot ® 


D ,a '^e v . a { v. \ 
E er "dr -F® r 3r\T®) 




ea« iy i dv e civ 0 dv 0 cote\ 
>-_ « r\sl ne d0 + de 3® / 


H 

I 


-rf ’cos 0 


( 2 . 1 ) 


-rf cosec ose 


( 2 . 2 ) 


rfsin®cos< 20 ) 


rfcos2©sin0 rfsln2esin20 



* 2 faln 6 cos 0 - 4 fcos©cos 20 


- 2 fs)ndco $0 ♦ 4 f cosecos 20 


-rf'cosecos® rf'sinecos20 


♦rf *cosesin0 rf'cos 2 esin 0 rf 'sln20sln20 


-fsin26sln0 2 f (l*cos*©)sin20 



fau 


y m 2( rf ' >*. 


d-sln’esin’0) 


8 f * ( 4 cos a e* 
2 (rf * ) *0j (0,0) sin 4 esln’20) 


ef’sin'e a 


♦ 8 (rf ’ )’sin a e 
2 « ( l-sl n*esi n*20) 


4 


Table 1.2 Eigenfunctions (l,m,n) numerically calculated for 
•140 radial zones in the solar convection zone. E DOT=viscous 
energy change=E v in eq. 1.3; K.E.=kinetic energyrSIGMA PRIME 
damping rate. 


E DOT K.E. SIGMA PRIME 

EIGEN 1,1,0 - 




2 . 501470E + 11 

7 . 532208E + 17 

1 . 235631E-05 

EIGEN 

1,1,1 - 

E DOT 
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SIGMA PRIME 



2 . 955946E + 14 
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3 . 480709E-05 
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modulations at the rotation rates, repetition rates and beat 
frequencies of envelope r-modes (in the CZ). To first order, the 
® r-modes- oscil late and rotate at rates independent of all solar 

interor conditions except for the rotation of the fluid in which 
they are embedded. Thus the r-mode rotation rates are as predictible 
as the g-modes. Each r-mode sequence is determined by a single 
free parameter: a mean rotation rate of the cavity in which the 
# modes are trapped. This implies regular but complex behavior in 

solar characteristics modulated by the r-modes. In summary, some 
periods of r-modes have been tabulated. Viscous damping is strong 
and removes one-third or more of the mode energy each oscillation 
period. 

** Table 1.3. Sidereal rotation of r-modes in the convection zone 

and the repetition of identical flow fields as viewed from Earth 
(Wolff & Blizard, 1986). 

1 m ROTATION REPETITION, l> lm 

® Rate Period Rate Period 
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2.1 FAST FOURIER TRANSFORMS OF THE ZURICH DAILY SUNSPOT NUMBER 
AND GREENWICH DAILY WHOLE SPOT AREA 

A -fast Fourier transform (FFT) of daily Zurich sunspot 
number was made for eight sets of data, taken near peaks of solar 
cycles 14-21. The FFT accepts 2048 bits of data, so 2048 days 
were included, starting with 1 Jan except for the most recent 
cycle, when the data set began 16 Nov 1978, when data was 
available from the ACRIM irradiance satellite experiment. 

Thus the eight files contained 2048 days each, 
with a starting date of: 


File 1 16 Nov 1978 
File 2 01 Jan 1966 
File 3 01 Jan 1956 
File 4 01 Jan 1945 


File 5 01 Jan 1935 
File 6 01 Jan 1925 
File 7 01 Jan 1915 
File 8 01 Jan 1904 


For a sequence of 2048 days, the FFT will display periods 
between 2 days and 1024 days (2.8 yrs) due to the pairing of the 
data. Prominent in most of the spectra is the 27-day peak. 

Periods longer than 160 days would not be resolved. 

The Greenwich Photoheliographic Results, 1874-1981, (GPR) , 
were obtained on magnetic tape from the National Geophysical Data 
Center of NOAA. In addition, an updated version of GPR 1982-1984 
was obtained from Douglas Hoyt of NCAR. The GPR tape was loaded 
onto a VAX and pertinent columns extracted (date, whole spot 
area). A program was written to average data over four days. 
Another program was written to fill in zeros of whole spot area 
(when no observations were made), by averaging the values of the 
day before and day after. Because of large data gaps in early 
years, the period 1904-1984 was selected for study with fast 
Fourier transform (FFT) (Fig. 2.0). 

FFTs were run on eight data sets with the same starting dates 
as for the Zurich sunspot number. In this case, the plots were on 
a single frame, rather than four frames for the previous data 
sets, since there were 512 points rather than 2048 (4-day 
averages) (Figs. 2. 1-2. 4). The mean of the eight spectra showed 
very little detail. For the 4-day averages, the limits on the FFT 
are 8 days and 5.6 years, but in practise the lower limit is 16 
days because of the mirror image form of the Graf plot. Because 
the four-day smoothing resulted in improved spectra, it was 
applied to the Zurich sunspot number data, which were replotted 
(Figs. 2. 5-2. 8) . 
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Power Spectrum, Daily Whole 
Spot Area: 4-day averages 
Cycle 16: Jan 1925-Aug 1930 
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b 

Power Spectrum, Daily Whole 
Spot Area: 4-day averages 
Cycle 17: Jan 1935-Aug 1940 



Period, days 


Fig. 2.2 Fast Fourier transforms of Greenwich projected whole 
spot area (4-day averages): solar cycles 16 and 17 (a, b). 
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Fig. 2.3 Fast Fourier transforms of Greenwich projected whole 
spot area (4-day averages): solar cycles 18 and 19(a,b). 
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Fie. 2.4 Fast Fourier transforms of Greenwich projected whole 
spot area (4-day averages): solar cycles 20 and 21 (a, b). 




2.5. Fast Fourier transforms of Zurich daily sunspot number 
(4-day averages) for solar cycle 14 (a) and 15 (b). 
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Fig. 2.6. Fast Fourier transforms of Zurich daily sunspot number 
(4-day averages) for solar cycle 16 (a) and 17 (b) . 
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Power Spectrum, Zurich 
Daily Sunspot Ho. (4-day ave) 
Cycle 18: Jan 1945-Aug 1950 
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Fig. 2.7. Fast Fourier transfo 
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Power Spectrum, Zurich 
Daily Sunspot No. (4-day ave) 
Cycle 21: Nov 1,978-Jul 1984 
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Fie. 2.8. Fast Fourier transforms of Zurich daily sunspot number 
(4-day averages) for solar cycle 20 (a> and 2l (b) . 
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3 EVALUATION OF FFT'S WITH WOLFF & HICKEY (1987) MODEL 

Each theoretical frequency having measurable power which 
stands out prominently in any spectra is listed in Table 3.1. 
Frequencies prominent in the earlier cycles were primarily caused 
by r-modes. Only one of the 7 cases where a line was clearly more 
prominent in the earlier cycles did not involve r-modes; 5 of the 
7 were exclusively caused by r-modes. A different situation was 
found in the last cycle (21) (later in the solar cycle), where 
only 2 of the 8 cases did not involve g-modes. This is consistent 
with decaying r-modes which would be expected in the Sun. Wolff 
and Blizard (1986) found that turbulent viscosity would damp r- 
modes in a few months, and postulated that they might be excited 
irregularly by large convection events. Occasional injections of 
fresh energy would prolong r-mode oscillations. In any case, the 
above measurement showed r-modes to be a stronger influence in 
the first few years of solar cycles rather than later in the 
cycle, when their lessened influence allowed the presumably more 
stable g-modes to have relative prominence (see Fig. 2.0 & 3.1). 

Figure 3.1 shows a Fourier spectrum of projected sunspot 
area (four-day averages). Below it lies the model lines from 
Wolff & Hickey (1987) and Table 3.1. The region between zero and 
50 nHz contains many theoretical lines due to high-1 modes, 
indicated here only as a group (labeled H). The data is 
unresolved but the high power measured there is expected. At 
frequencies above 50 nHz, there is a remarkable similarity 
between theory and the observed sunspot area spectra. Predicted 
lines occur where there is much observed power and minima in 
power tend to lie where there are no theoretical lines. All four 
of the prominent power maxima near 310, 380, 450 and 530 nHz are 
associated with the lowest angular harmonic. The first three 
correspond to beats involving the 1=1 r-mode. The fourth involves 
the 1=1 g-mode and the feature H. This supports the expectation 
that the lowest modes are the most influential on variations of 
sunspot area (and other indicators of solar activity). 

Only three peaks (100, 170 and 295 nHz) occur near no 
predicted line. Finally, considerable power is seen in the 7 
bands which all arise from feature H. All seven lie where there 
is considerable observed power and 5 of the 7 include observed 
maxima. If the model is correct, it implies that the combined 
effect of all high harmonic r-modes can be comparable to one low 
harmonic in modulating sunspot area. Overall, the agreement seen 
in Table 3.1 and Figure 3.1 is much better than one could expect 
by chance, since there were effectively fewer than two free 
parameters to adjust the 28 model frequencies (Wolff & Hickey, 
1987) . 
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| 4.1 MAGNITUDE OF EXTERNAL TORQUE 

We wish to calculate if the horizontal component of the 
force exerted by a planetary torque, does work on the photo- 
spheric layers at a rate comparable to 10“® L Q (the mean rate of 
energy seen released by solar activity). Assume that F h is 
constant and acts for about 5 days (1/6 rotation of the sun) on a 
* volume element of the photosphere, accelerating it away from the 

equator with no resisting force (for simpl icity) (Fig. 4.1). Then 
the horizontal velocity is v^ * (F^)t, where F^ is in units of 
force/mass, then (kinetic energy/ volume) = 1/2 pv^ . In choosing 
9 , we do not go into the convecting interior because these layers 
I may have a stiffness that resists the torque-induced flow, since 

they experience strong horizontal diffusion of momentum. This 
would give an upper limit to an energy source fed into the photo- 
sphere. The long-term mean rate of energy deposition would be: 

(vol .equatorial bulge in photosphere) xp/ 2 [v^ (for t=5 days)] * 

9 15 days 

The force of an external body on the sun may be calculated from 
McCul lagh's formula, which expresses the gravitaional potential 
of any body at a distance R and with a mass M(see section 3.1): 

I V = GM/R + GM (A+B+C-3I) /2R 3 + higher order terms (4.1.1) 

where A, B and C are the principal moments of inertia of the 
Sun. It is usual to write J 2 =3/2 (A-C) /MR Q 2 , where J 2 is the 
quadrupole moment. Since we assume that the sun is axisymmetric, 

I to a first approximation, A=B, and (A-C)/C is approximately 10 “ 5 

(Hill, 1975). The torque of an external body is obtained by 
finding the force from the above equation (4.1.1) by differen- 
tiation; since the first term vanishes, the torque is then FxR, 
or, (see Fig. 4,2 for sin 9), 

V = 3GM (A-C) sin e/R 3 (4.1.2) 

^ Using I =3xl0^kg-m , and G=6. 67x10” m 2/kg , 

I‘ = 36.67xl0 _11 Nm 2 /kg 2 ) (3xl0 4b kg-m 2 ) (10" 5 )MjSine/R 3 (4.1.3) 

Substituting for Jupiter, M j=19.1xl0 27 kg, Rj=7.78xl0 11 m, and 
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sin 6°=0.104, then the torque due to Jupiter is =24.0xl0 24 
(Joules). For other planets, p =20 . 0xl0^ 6 n-m 2 /kg^ (MpSin©/R 3 
-• Next, calculate the volume V b of the oblate equatorial 
of the Sun: 


N-m 
) . 

belt 


where AR/R Q 


V b = (2trR e )4R4H (4.1.4) 

IQ" 5 and, arbitrarily, AH = 2xl0®m, then 


= 6.2xl0 2 *m^ 


(4.1.5) 


Now, find the mass of the oblate belt, using a density for the 
photosphere p = 3xl0“ 5 kg/m 3 (Abell, 3rd Ed). Then, 

M b = P V b = 1 8.6xl0 16 kg (4.1.6) 

Consider the torque from one planet, say, Jupiter (equation 
4.1.3 and so on). Find the average horizontal force from the 
quadrupole term, 

F h =Fj/R 0 = (24.0xl0 24 N-m) /7 . 0xl0 8 m=3 . 42xl0 16 N (4.1.7) 

Calculate the horizontal velocity v b =(F b )t, where t=5 days, and 
F b is the force per unit mass (divide by mass of belt, M b ) : 

v b = F h /M b (5 days) = 79.4 km/s = 7.94xl0 4 m/s (4.1.8) 

which is the same order of magnitude as the equatorial velocity 
of the sun, equal to 2 km/sec. 

Now calculate the kinetic energy per unit time, 

(KE/sec) = V b x(l/2)p(v h ) 2 /15 days 

= 45.2xl0 19 N-m/s = 4 5. 2xl0 19 J/s (4.1.9) 

Compare this figure to the energy of solar activity, 10“^L Q , 
which is 4xl0 2 ®j/sec (Allen, 1973). Thus, 

(KE/sec) = 45.2xl0 19 J/s > (lO -6 ) (4xl0 20 J/s (4.1.10) 

(KE/sec) > 10" 6 L q 

So the energy input is much greater than the estimated energy of 
solar activity. 



4.2 TINE SPAN OP TORQUE IN EXTERNAL SECTOR OF SUN 

For, the case of two planets, the problem is to find how long 
both planets are located in a sector of given angular width as 
viewed from the sun, near heliographic conjunction. Analytically, 
one uses angular velocity, w, or the equivalent, 2II/P, where P is 
the planet period in convenient units, such as days. Then start 
from the two innermost planets and proceed outwards. Let t be the 
time both Mercury and Venus are located in a given sector of 
heliocentric longitude, 6, then, 

46 = (Aw)t so t = 46/4w = 46/ (w 1 - w 2 ) (4.2.1) 

or using periods, 

t = Ae/2H(P 1 P 2 ) / (P 2 -P 1 ) (4.2.2) 

when angles are expressed in radians. If the angles are in 
degrees, multiply by 360°/2H. If one uses an angular sector of 
1/2 radian (28.65°), then the time in equation (4.3.2) becomes 

t = 1 (87.969 d) (224.70 d)_ 5 11.50 days (4.2.3) 

4n (224.701 - 87 . 9 69 ) d 

Similarly for all other combinations of two planets. Using 
sidereal periods in days (Allen, 1973) one obtains the values in 
Table 4.1 A graph can be used to check results. Find the 
degrees/day each planet moves in longitude (Astron. Almanac) and 
plot on a solar system chart. Diagram the time taken for a given 
planet to move 28.65° in heliocentric longitude, and compare with 
another planet when they are in conjunction (see Fig. 4.2.1, 
Tables 4.1 & 4.2). 

The axial rotation of the sun does not enter here. Using the 
axial angular velocity of the sun, w Q , and the two planetary 
revolutionary velocities, w-^ &w 2 , it can be shown that w Q cancels 
out (where w^'and w 2 ' are synodic velocities): 

4,w = Wj_* - w 2 ' but Wj' * w 0 " w i 

so, 

4w = (w G - w^ - (w 0 - w 2 ) = w x - w 2 (4.2.4) 


Table 4.3 shows the "synodic” periods of the planets as viewed 
from the subsolar point on the sun's surface. 

Finally, Table 4.4 lists the relative amplitude of 
precessional torque of the various planets. 
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Fig. 4 . 2. 1. Plan view of solar system, showing planet positions bv 1976 date 
Time span when 2 planets are within a sector of half radian near conjunction. 
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Table 4.1 Days that two given planets are within a sector of 4 radian 
viewed from the sun. (Allen. 1973). 


- 

P(days)* 

Venus Earth Mars Jupiter Saturn Uranus Neptune 

Her 

87.969 

11.504 9.221 8.028 7.145 7.06 7.02 7.01 

ven 

224.701 

46.47 26.572 18.86 18.26 18.01 17.9 

Ear 

365.256 

62.065 31.74 30.09 29.42 29.2 

Mar 

686.780 

64,95 58.40 55.92 55.3 

Jup 

4332.589 

577.21 401.46 371.5 

Sat 

10759.22 

1318.50 1042.6 

Ura 

30685.4 

4834 . 1 

Nep 

60189. 



it sidereal 

Table 4.2 Days that two planets are within 4 radian, heliocentric 
longitude, viewed from the sun; four solutions (Astron Almanac) 


Conj 

1976 dates 

Planets 

Calc. 

days 

Graph P * /41T 

Daily M 

deg/day 

20 Jul 

Me 

V 

11.5 

11.0 

11.5 

7.78 

5.28-1.6° 

22 Sep 

Me 

E 

9.22 

11.5 

9.21 

8.32 

4.44-1.0° 

05 Jul 

Me 

J 

7.15 

9.0 

7.15 

4.B3 

6.06-0.09° 

15 Jun« 

V 

E 

46.5 

44.0 

46.5 

46.2 

1.62-1.0° 

19 May 

V 

J 

* 18.9 

20.0 

22.4 

18.9 

1.60-0.09° 

18 Nov 

E 

J 

31.7 

30.5 

31.8 

31.5 

1.00-0.09° 


■superior conjunction 


Table 4.3. Periods of the planets, viewed from the subsolar 
point. The sidereal rotation period of the sun is taken as 
25.80 days (Wolff. 1976). (Allen. 1973) 



Sidereal Period 

Svnodic Period 

Mercury 

87.969 d 

36.51 d 

Venus 

224.701 

29.15 

Earth 

365.256 

27.76 

Mars 

686.790 

26.81 

Jupiter 

4332.589 

25.95 

Saturn 

10759.22 

25.86 

Uranus 

30685.4 

25.82 

Neptune 

60189. 

25.81 


Table 4.4 Relative amplitude 
the sun. (Alle n. 1973). J - 


PLANET 

G m (10 mVsec 2 

Mercury 

22031.8 

Venus 

324858.6 

Earth +Moon 

403503.22 

Mars 

42828.44 

Jupiter 

126686.9 x 10 3 

Saturn 

37938.5 x 10 3 

Uranus 

585 x 10 4 

Neptune 

686 x 10 4 


of precessional torque of planets on 
-4 x 10' s (Hill & Stebbins, 1975). 


T 


G m/r* 


G m sin 28/r 3 


1.3 

X 

10*'’ 

1.5 x 10~ 

2.6 

X 

lO"'* 1 

3.4 x 

1.2 

X 

io-* 

3.0 x 

3.7 

X 

IQ' 1 ' 

0 . 072 

2.7 

X 

10'” 

5.6 X 10" 

1.3 

X 

10'*' 

0.24 

2.5 

X 

10 

0.0055 " 

7.5 

X 

io' 15 

0.00016 * 


.-20 


.-20 
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Fig.4.3-lfa> Two point masses are situated on the y-axis at ia. 
Find the gravitational field along the x-axis, when x>>a. 

(b) The half-value g^ of 0 max - <See text). 

4.3 THE COMBINED GRAVITY FIELD OF TWO PLANETS ACTING ON THE SDN 

We wish to calculate the angular separation in longitude of 
two gravitational bodies/ such that the gravitational field at a 
large distance perpendicular to the line joining them is one-half 
of the value of the two bodies superimposed (in line). This would 
approximate two planets in conjunction, then separated by an 
angle 0 (both have the same mass and are at the same distance as 
a first approximation) (See Fig. 4.3.1). 

To find the half-value of g x , set g x = (l/2)g max . Take a 
ratio _ _ , 

q„ = 1 = -Gmx/ (x 2 + a 2 ) 3 ' 2 (4.3.1) 

W 2 -Gmx/(x 2 ) 3/2 

Solve for x, calculate tan 0 = a/x = 0.7665; 0 = 37.47°. 

Thus two planets of the same mass and in the same orbit separated 
by 37.47° longitude have one-half of the gravitational force that 
they would have if at the same longitude (in conjunction) . 

In actuality, the duration of effect of two planets in 
conjunction on sunspot number extends to *_ 5 days, in the case of 
Venus and the Earth (Link et al, 1954), or about 10°, not 37°, so 
the Gaussian or cosine formula to be used should be a much 
steeper function. 
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4.4 NORTH-SOOTH ASYMMETRY OF SOLAR ACTIVITY 

Solar activity is not equally divided between the northern 
and southern hemispheres of the sun. Before 1850, sunspots were 
more numerous in the northern hemisphere. From 1865 to 1910 they 
predominated in the southern hemisphere. From 1910 to 1969 there 
was a northern excess, with a strong asymmetry in cycle 20 
(Dodson & Hedeman, 1970) . Apparently in the 1980's the northern 
excess may have come to an end. 

If sunspots (and associated solar activity) are influenced 
by planetary positions, then their slow swing in N/S latitude may 
be related to the similar slow latitude swing in the combined 
influence of al 1 planets which exert an effect on the sun. Any 
such influence should refer to planetary latitude with respect to 
the solar equator. The slowest important planets, Uranus and 
Neptune, reached their greatest latitude north only 10 years 
apart in 1672 and 1682, respectively, in the solar equatorial 
system. The most extreme case of N/S sunspot asymmetry occurred 
between 1672 and 1704, when no northern spots were observed. 

During this entire time period two or more outer planets were 
north of the solar equator. Two years before Uranus reached its 
descending node, Saturn passed its ascending node. 

The orbital period of Uranus is 84 years. The conjunction 
period of Saturn and Neptune is 36 years, or about half the 
observed long N/S period of sunspots. The Neptune orbital period 
is 164 years, or about twice the observed N/S sunspot period (See 
Figs. 4.4.1 & 4.4.2 and Table 4.4.1). 

Various observers find a long period N/S oscillation of 
sunspot positions of 55 years (Berdischevskaya, 1967) of 70 and 
84 years (Brunner-Hagger & Liepert, 1944), and 70 years (Wolbach, 
1960) (White & Trotter, 1971), depending upon how it is measured. 

Shown in Table 4.4.1 is the latitude of three outer planets by 
decades for the entire historical record of sunspot measurement. 

In Fig. 4.4.2 is shown a comparison of the N/S asymmetry of 
sunspots for the period 1943-1977 (Fracastoro & Marocchi, 1978), 
and the latitude of four outer planets relative to the solar 
equator. 

Also relevant is the latitudinal position of the sun 
center (CS) relative to the Center of Mass of the Solar System 
(Petrova et al, 1978). The outer planets contribute the following 
percentages to the displacement of the sun: Jupiter: 49%, Saturn: 
27%, Uranus: 8%, Neptune: 15%, and Pluto: 1%. The vector quantity 
(CS-CMSS) is related to the sun surface and sun center Coriolis 
acceleration (Pimm & Bjorn, 1969). 

The K-corona and the zodiacal light may rotate at an 
inclination to the solar equator. Evidence is seen from 
observations of the K-corona and zodiacal light, especially when 
observations are made near to the axis of the nodes of the 
planetary system on the solar equator. The inner K-corona is 
inclined about 3° to the solar equatorial plane. Farther out the 
zodiacal cloud is inclined between 4° and 5.7° (Misconi, 1977). 
There would be a resulting non-uniform shear around the solar 
equator, which might engender oscillations. Such oscillations may 
have been observed, with periods of half a year (180.5 days) and 
half of Jupiter's synodic period (199.3 days) (Czada, 1983). 


Fig. 4.4.1. Relative spotted area of the north and south hemispheres , left, 
expressed by the ratio q* (n-s)/(n*s) , compared to latitude of outer planets 
above or below the solar equator (Wolbach, 1962). 



it l ?0th century year"** | 

« US U1 U9 51 53 55 57 59 61 63 65 67 69 71 73 75 77 

Ficr. 4,4.2 Monthly values of parameter &. the asymmetry between 
the activity in the north and south hemispheres 1943-1977 (left) 
(Fracastoro a Marocchi, 1978), compared to latitude of outer 
planets relative to the solar equator (right). 

Table 4.41Longl tude of the outer planets, Jan l or year given, 
relative position north or south of the solar equator, and 
approximate Latitude <B) from the solar equator. 


Jan l 
Year 


1840 

1850 

1860 

1870 

1880 

1890 

1900 

1910 

1920 

1930 

1940 

1930 

1960 

1970 

1980 

1990 

2000 


353 * *6 
128 * - 3 * 

216* -3* 

3 * * 6 * 

135 * - 4 * 

240 * - 0 . 3 * 
350 * * 6 * 

130 * - 4 * 

230 * - 0 . 5 * 

3* *6* 

140 * - 4 . 3 * 
254 * * 0 . 5 * 
10* *1* 

150* -5* 

270 * * 1 * 

25 * ♦»* 

160 * - 6 * 

280* *1" S 
40* *1* S 
179* -5* S 
286* *1.5* N 
47* *0.2* N 


Neptune 
N/S X 8 


260* *1 
303 * + 3 * 
343 * * 6 * 
26 * * 3 * 
70 * * 0.5 


-5* 
-3* 
.3 

286 * * 1 “ 
330* *3* 


135 * - 4 * 
185 * - 6 * 
230 * - 1 * 
271 * -l* 
313 * * 4 * 



N 

247 * 

0 * 

N 

269 * 

♦ 1 * 

N 

291 * 

♦ 2 * 

N 

310 * 

♦ 4 . 5 * 

N 

335 * 

♦ 6 * 

N 

0 * 

* 5 * 

N 

19 * 

♦ 3 . 3 * 

N 

45 * 

* 1 . 5 * 

' 5 

64 * 

- 0 . 2 * 

S 

90 * 

- 1 * 

S 

108 * 

- 2 . 5 * 

S 

135 * 

- 5 * 

s 

152 * 

- 6 * 

s 

175 * 

- 3 . 5 * 

s 

195 * 

- 4 * 

s 

215 * 

- 2 * 

s 

238 * 

- 0 . 2 * 

N 

256 * 

♦ 0 . 5 * 

N 

283 * 

♦ 2 * 

N 

305 * 

♦ 4 * 
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